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We demonstrate that the 'microscopic' field theory representation, directly derived from the cor- 
responding master equation, fails to adequately capture the continuous nonequilibrium phase tran- 
sition of the Pair Contact Process with Diffusion (PCPD). The ensuing renormalization group (RG) 
flow equations do not allow for a stable fixed point in the parameter region that is accessible by 
the physical initial conditions. There exists a stable RG fixed point outside this regime, but the 
resulting scaling exponents, in conjunction with the predicted particle anticorrelations at the critical 
point, would be in contradiction with the positivity of the equal-time mean-square particle number 
fluctuations. We conclude that a more coarse-grained effective held theory approach is required to 
elucidate the critical properties of the PCPD. 



PACS numbers: 05.40.-a, 64.60.Ak, 64.60.Ht, 82.20.- 



I. INTRODUCTION 

Phase transitions between different nonequilibrium 
steady states are frequently encountered in nature, and 
determining the associated critical properties is an im- 
portant issue. Unfortunately, compared with the sit- 
uation in thermal equilibrium, a full classification of 
nonequilibrium phase transitions is still in its infancy. 
We shall focus here on a particular subclass of nonequilib- 
rium transitions which separate an 'active' phase, char- 
acterized by a fluctuating order parameter 0(r, t) with 
nonzero average (0), from an absorbing state wherein 
((/>) = 0. In the thermodynamic limit, all degrees of free- 
dom remain strictly frozen in such an inactive, absorbing 
phase Pj. 

The universality classes of such transitions are con- 
veniently studied in the framework of reaction-diffusion 
processes, even though other descriptions abound (sur- 
face growth, self-organized criticality) 2j. The most 
prominent representative of absorbing state transitions 
is the contact process (CP), for under quite generic con- 
ditions, namely spatially and temporally local micro- 
scopic dynamics, and the absence of coupling to other 
slow fields (thus excluding quenched disorder and the 
presence of conservation laws) , active to absorbing state 
transitions fall into the CP universality class with scal- 
ing exponents that also describe critical directed perco- 
lation (DP) clusters 0, fl- Yet, the very fact that de- 
spite considerable effort hardly any experiments have to 
date unambiguously identified the CP/DP critical ex- 
ponents, hints at the prevalence of other universality 
classes. In simulations, the ^ parity- conserving'' (PC) uni- 
versality class is also prominent: represented by one- 
dimensional branching and annihilating random walks 
(BAW) A ^ {m -\- 1)A, A + A ^ with even m, it 



is characterized by local particle number parity conserva- 
tion. In contrast, the phase transition in low-dimensional 
BAW with odd m is governed by DP exponents 0, ■ 

Novel critical behavior is to be expected when all in- 
volved reactions require the presence of neighboring par- 
ticle pairs 01 . The pair contact process with diffusion 
(PCPD, or annihilation/fission model |^J) is conveniently 
defined through the microscopic reaction rules 



A + A 



A + A'^A, A + A^A + A + A (1) 



(the presence of either pair annihilation ~ /i or coag- 
ulation ~ /i' suffices), supplemented with particle hop- 
ping (diffusion constant D) subject to mutual exclusion. 
The latter is crucial for the existence of a well-defined 
active phase and continuous transition. For without 
restrictions on the occupation number per lattice site, 
the particle density diverges within a finite time when 
(T > (Tc = -|- /Lt' 1^. In the inactive phase, however, 
site exclusion should not be relevant. It is then easily 
seen that the absorbing state of the PCPD (as in the PC 
universality class (5|) is governed by the algebraic den- 
sity decay of the pure pair annihilation process 0, viz. 
{(j){t)) - in dimensions d > 2, {(j){t)) ~ t''^/'^ for 
d < 2, and {4>{t)) ^ Int at the (upper) critical dimen- 
sion dc = 2. In contrast, in the CP/DP universality class, 
the inactive phase is characterized by exponential particle 
decay and correlations. Recall that here the branching 
processes merely require the presence of a single particle: 
the third reaction in would simply be replaced with 
A ^ A + A. Site exclusion is not crucial in this case, 
as long as the pair annihilation or coagulation reactions 
are included. Alternatively, the combined first-order re- 
actions yl — > and A ^ A -\- A with site exclusion yield 
a CP/DP continuous phase transition as well. 
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Holding the rates fj, and /i' fixed, there is a critical pro- 
duction rate ac at which the transition between the active 
nonequilibrium steady state and the absorbing phase oc- 
curs. It is a central issue, in an effort to classify nonequi- 
librium phase transitions, to clarify the precise manner 
in which the particle production mechanism defines the 
properties of both the absorbing state and the univer- 
sality class of the transition, i.e., how it affects scaling 
properties in the vicinity of the critical point. Numeri- 
cal investigations of the PCPD started with Ref. |8|. It 
almost constitutes an euphemism to state that this and 
the subsequent flurry of numerical work |j^-|l9|| have re- 
vealed conflicting views (see Ref. |^ for a comprehensive 
overview): For not only are the precise numerical val- 
ues of the critical exponents still being debated to this 
day, but even more striking, the very issue of the PCPD 
universality class has remained controversial. Essentially 
three scenarios have been put forward: Either the tran- 
sition defines a novel independent universality class that 
is yet to be characterized, or it belongs to the CP/DP, 
or even to the PC class (the latter perhaps becoming less 
likely with improving simulation accuracy). In addition, 
the emergence of different scaling properties depending 
on the value of the diffusion rate has been claimed. 

This inconclusive numerical situation clearly calls for 
analytical approaches to provide further understanding 
of the elusive continuous phase transition of the PCPD 
with restricted particle occupation numbers. A natural 
starting point is the standard field-theoretic representa- 
tion of reaction-diffusion systems that can be derived di- 
rectly from the corresponding classical master equation 
[21I. Specifically, dynamical renormalization group (RG) 
studies based on such 'microscopic' field theories were, 
e.g., successfully employed to diffusion-limited annihila- 
tion and even-offspring BAW |^, as well as the inac- 
tive phase of the PCPD without site occupation restric- 
tions 1^. In either case, particle anizcorrelations gov- 
ern the asymptotic scaling regime, as opposed to the 
typical clustering behavior in the CP/DP universality 
class (that includes odd-offspring BAW) (3|. Here an- 
other coarse-graining step takes the original 'microscopic' 
master equation representation to Reggeon field theory, 
equivalent to a Langevin equation with 'square-root' mul- 
tiplicative noise, that serves as the appropriate effective 
action for the CP/DP critical properties. Thus, one 
would hope that the continuous nonequilibrium phase 
transition in the PCPD with restricted site occupations 
should be amenable to these powerful tools as well. 

Yet it was only recently demonstrated how site ex- 
clusions can be consistently incorporated into the mas- 
ter equation field theory formalism 22]. This paper re- 
ports our study of the ensuing action for the PCPD, 
constructed in Sec. ^1 carefully taking into account the 
higher-order reactions that become generated through 
fluctuations, i.e., successive particle production processes 
m. We shall derive and discuss the ensuing RG flow 
equations in Sec. IIIII We will demonstrate that there 
exists in fact no stable RG flxed point in the physically 



accessible parameter space of the model (wherein all reac- 
tion rates are non-negative). Remarkably therefore, the 
microscopic field theory, albeit directly derived from the 
master equation, is not capable of capturing the PCPD 
critical point. We interpret the appearance of runaway 
flows as an indication that a crucial ingredient was obvi- 
ously left out when the (naive) continuum limit was taken 
j23l |. An appropriate effective coarse-grained description 
might require the explicit introduction of separate den- 
sity fields for the 'inert' solitary random walkers and the 
clustered particles, respectively, akin to the explicit treat- 
ment of particle pairs and singlets in Ref. ^10]. The RG 
flow equations do however allow for a fixed point outside 
the physical region. In Sec.|^we compute the associated 
critical exponents to second order in an expansion around 
the upper critical dimension dc — 2, and moreover estab- 
lish exact scaling relations. But we shall see that the ac- 
tual exponent values, when combined with the predicted 
particle anticorrelations at the critical point, violate the 
positivity of the equal-time mean-square particle fluctu- 
ations. Hence this RG fixed point is clearly unphysical. 
In conclusion, the construction of a consistent field the- 
ory description of the PCPD remains an open problem. 
Similar issues arise also for the closely related models 
that involve solely particle triplet or quadruplet reactions 



II. MASTER EQUATION FIELD THEORY 
REPRESENTATION OF THE PCPD 



The classical master equation kinetics of particles sub- 
ject to diffusion and local 'chemical' reactions can be 
mapped onto a field theory action following standard 
procedures However, for the density to remain 

bounded in the processes at arbitrary values of the 
reaction rates, specifically in the active phase, it is neces- 
sary to introduce a growth- limiting process. In most nu- 
merical simulations this is achieved by further imposing 
mutual exclusion between particles. Analytical progress 
therefore requires a consistent incorporation of the ex- 
clusion constraint. To this end, we follow Ref. [U, and 
write down the resulting action corresponding to the pro- 
cesses (QJ on a (for the sake of notational simplicity one- 
dimensional) lattice (sites i): 



S[{$i, (j>i}] Jdt dt<l>i + {fi(<j>i ~ 1 



' 2 
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X (picpi+i > (Pi (Pi+i e 



The time-dependent fields 4ii{t) and 4)i{t) here originate 
from a coherent-state representation employing hosonic 
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creation and annihilation operators |3i|211. The exclusion 
constraints are encoded in the exponential terms [2^ . 
and the unrestricted model is recovered when all these 
exponentials are replaced with unity. 

Thus far, the action (0) constitutes an exact repre- 
sentation of the microscopic processes (0. In order to 
proceed to a continuum field theory, which should suf- 
fice to describe the large-scale, long-time behavior in the 
vicinity of a critical point, we add a diffusion term (for 
which we ignore the site occupation restrictions '26^1) and 
take the (naive) continuum limit (now in d dimensions) 
<j>i{t) — > <j>{r,t) and <j>i{t) a'^(l){r,t), with a denoting 
the original lattice spacing, such that 0(r, t) remains di- 
mensionless. This yields 

5'[0,0] = Jd'^x Jdt l^4>{dt- DV^)4> 



^2 



(3) 



with a microscopic inverse density scale u ~ a"*. 

The corresponding classical rate equation (augmented 
with diffusion) is readily obtained by solving for the sta- 
tionarity conditions (55[0, 0]/(5(/) — 0, which, as a conse- 
quence of probability conservation in the master equa- 
tion, is always satisfied by </> = I, and 6S[(I>, (f)]/5(/) = 0, 
which then results in 

{dt -DV^) 0(r, t) = _ (2^ + ^] g-^- * . 

(4) 

In contrast with the unrestricted model (where w = 0), 
this mean-field equation allows for an active state with 
a finite particle density (j) = ln(cr/crc), provided a > 
CTc = 2/i + /x'. Near the now well-defined critical point at 
CTr, we obtain 



0(a) 



^[(a/a,)-l]^(a-a,)^ P=l 



(5) 



In the absorbing phase (cr < CTc), the site restrictions 
do not matter, and the density decays algebraically as in 
pure annihilation or coagulation, (j){t) ~ [{ac — <T)t]~^ . At 
the critical point, this becomes replaced with the slower 
power law 



(bit) ~ (w,t)-i/2 _r 



<5= 1/2 



(6) 



This relation already shows that the scaling of the pa- 
rameter V determines the critical exponents. Moreover, 
scaling analysis tells us that the static correlation length 
diverges upon approaching the critical point from the 
active side according to ^(cr) ~ (cr — with = f , 
whereas the characteristic time scales as tc ^ with 
diffusive dynamic critical exponent z — 2. 

While we clearly need to retain the exclusion parame- 
ter in order to describe the continuous phase transition 
occuring at A oc cr — ctc = 0, we also note that an expan- 
sion to first order in v suffices. More technically, since 



the field (f){r, t) scales as a particle density, the scaling 
dimension of the exclusion parameter is [v] = k~'^, where 
K represents an arbitrary momentum scale. Superficially, 
therefore, v constitutes an irrelevant coupling that flows 
to zero under scale transformations. We may thus ex- 
pand the exponentials in the action keeping only the 
lowest-order contributions, which leads to additional in- 
teraction terms. Upon at last performing the field shift 
(/>(r, t) — l-|-0(r, t) (whereby final-time contributions, not 
explicitly listed here, become eliminated |27j|'). we arrive 
at an action of the form 



dx / dt 



+A(0) 03 



(7) 



where we have defined the generating functions 



G{x)=J29p^r A(x)=^Apa;f. (8) 

p>i p>i 

Note that probability conservation implies G(0) =0 and 
A(0) — 0. Microscopically, we identify gi = 2/i + /i' — cr 
Ctc - CT, .g2 = M + - 2cr, ^3 = -cr, Ai = (3cr - 4/i - 2^j.')v, 
A2 = (9cr - e^i - 4:fi')v, A3 = (9ct - 2/i - 2fi')v, and 
A4 = 3cr V. However, at a coarse-grained level, fluctua- 
tions generate the entire sequence of particle production 
reactions 2A {n + 2)A (n > 1). For instance, two sub- 
sequent branching processes 2A — > 3A immediately lead 
to 2 A — > 4A, and so forth This effectively extends 
the sums in the functions ^ to all integer p. Unlike 
in conventional situations, we thus have to deal with an 
infinite number of vertices. Lastly we remark that in- 
troducing third-order annihilation reactions of the form 
3A kA [k — 0,1,2) also produces the terms in the 
second line of Eq. (TJ: Allowing for the back reactions 
of the particle production processes is equivalent to 'soft' 
site exclusions. 

With the previously introduced scaling dimensions of 
the fields [0] = hi^ and [0] = k'', we find [gp] — k^^"^ 
for all couplings in G(x), which suggests, as is then con- 
firmed by a careful analysis of Feynman diagrams, that 
dc — 2 constitutes the upper critical dimension here. 
Since [Ap] = k^'-^~'^\ the coefficients in the function A{x), 
which originates from site exclusion, are irrelevant near 
dc = 2. Yet because at least Ai is required to control the 
particle density in the active phase and thereby maintain 
a continuous transition, it cannot simply be omitted from 
Q. Once again though we arrive at the conclusion that 
terms of higher order in v [i.e., contributions ~ 0^ or 
higher in the action) need not be retained. But despite 
the presence of apparently infinitely many marginal cou- 
plings, the field theory Q remains renormalizable. This 
is best seen by recalling that the choice of the scaling di- 
mensions for the fields is actually arbitrary as long as the 
product [00] = k'^. Our theory thus contains a redundant 
variable l2a l that we fix conveniently as follows: Upon 



introducing recaled fields s = n'^^'^cf) and s 
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we obtain [gp] = k^-p'^/'^ and [Ap] = K2-(p+i)d/2^ 
sequently, the critical control parameter constitutes a 
relevant perturbation (for d < 4), since [51] = k^^''/^, 
whereas [32] = [Ai] = k'^~'^ indicating that both g2 and 
Ai are marginal at rfc = 2. Indeed, this procedure is 
consistent with the critical behavior according to Eq. © 
in two dimensions, which requires the density to scale 
~ t~^^^ ~ K rather than ~ ^ which is valid in the 
absorbing phase. All other couplings now acquire nega- 
tive scaling dimensions, and therefore become irrelevant 
for the leading scaling behavior. This leaves us with the 
reduced action 

S[s, s] = [dx [dt \s {dt ~ D\/^) s + giss^ 



+92 s^s^ + Aiss^ 



(9) 




Thus, the appropriate effective field theory for the critical 
point in fact only contains three nonlinear vertices. 



III. RENORMALIZATION AND RG FLOW 

It is instructive to proceed with the renormalization 
program based on either the field theory lO with in- 
finitely many marginal couplings and the reduced action 
One immediately notices that, to all orders in the 
perturbation expansion, the propagators do not become 
renormalized. Hence there is neither field nor diffusion 
constant renormalization, which already implies that the 
dynamic exponent in these field theories inevitably re- 
mains z = 2 exactly, at variance with present simula- 
tion data. Next, for the action (0 we define renormal- 
ized parameters according to gpC^ — ZpQp^iD k'^ , and 
Ap Ce = Zx^ ApR D with = (4^)~'^/2 r(l + e/2). 

The renormalization constants Zp and Z\ are deter- 
mined by means of dimensional regularization and min- 
imal subtraction by the condition that they absorb just 
the ultraviolet divergences appearing as poles in e = 2—d. 
After computing the RG (3 functions (3p = k df^gp^ (eval- 
uated in the unrenormalized theory) and upon introduc- 
ing the flow parameter £ = — ln(/ta) (where k is the run- 
ning momentum scale), we subsequently obtain the cor- 
responding RG flow equations for the running couplings, 
digpB.{i) = —I3p(t), and similarly for the ApR. 

We start with renormalizing the action J^. The ver- 
tex ^ Ai does not enter any Feynman diagrams that con- 
tribute to the renormalization of gi — Z\ gx^ D k}^^!"^ 
and 52 C'e = 'Z'2 52R -D k^. For the latter, we are therefore 
left with precisely the structure of the pure annihilation 
/ coagulation field theory 0, namely a geometric series 
of one-loop graphs, and hence arrive at the exact result 



Z-, 



Z-, 



I 



2.g2R 



(10) 



This leaves Zx^ in Ai = Zx^ Air D as the sole renor- 
malization constant to be actually determined anew here. 



FIG. 1: Qualitative flow diagram of the RG trajectories in 
the (gi, (;2)-plane. I denotes the stable fixed line; II and III 
represent the parts of the separatrix that might collapse in the 
limiting case with the negative half-axes. The hatched area 
indicates the absolutely unstable region. The dash-dotted line 
corresponds to the initial values of the model, with 2g\ — 32 = 
const. > 0. 



To two- loop order, we find [2 



6g2R. ^ 12giR 



2 3 



(11) 



From Eqs. H1U|) we obtain the RG flow equations 

di gi/2R(^) = [e - 2 52r(^)] ffi/2R(^) • (12) 

It follows that the sign of giR is invariant under the RG 
flow, and the critical giR — remains fixed. If both ^ir 
and (72R are negative, the flow according to Eqs. H12|) 
leads both running couplings towards —00 ('Wilson's 
gully'). Only if either parameter is initially positive, can 
the stable fixed line with arbitrary t/iR and (72R = e/2 
be reached. Thus, a separatrix must exist between the 
fully unstable region in parameter space and the basin of 
attraction of the fixed line. A qualitative sketch of the 
ensuing RG trajectories is depicted in Fig. ^ 

Before we proceed further with a discussion of the RG 
flow trajectories, let us consider the field theory j?)), and 
for the moment omit the irrelevant parameters Ap. One 
is then concerned with controlling the infinite number of 
marginally relevant vertices gp. This is most elegantly 
achieved by introducing the renormalized counterpart of 
the generating function Gr{x) = J2p>i 9pT^^^ 113 ■ 
anticipated, a careful analysis of the appropriate Feyn- 
man diagrams indeed shows that not only are the renor- 
malizations of the gp interwined, but also such couplings 
of arbitrarily high order become generated. For example, 
an explicit two-loop calculation results in 



2g2C^K " GgigsCeK 



De 



Dg2 e 



2Aglg4C^K 



-2e 



5 



(13) 



To one-loop order, elementary combinatorics yields Q 
Z,' = 1 - i E + 1) 3^^^^ 3^-^^^^ , (14) 

whence /3p = -e 5pR + X;j=i j (i + 1) 5i+iR ffp-i+iR- The 
ensuing infinite hierarchy of ome-loop flow equations for 
the (7pR is then efficiently recast into a single functional 
RG differential equation for GB.ix,£) [soj : 

d, Gn{x, e)=[e- dl Gr,(x, £)] Gr(x, i) . (15) 

Although we shall not explicitly make use of it, one may 
derive in a similar fashion the one-loop functional RG 
fiow equation for the generating function Ar that incor- 
porates all the couplings induced by the site exclusions: 

a,AR(x,^) = -2(l-e)AR(a;,£) (16) 
-GR(a;,£) 52Ar(x,£)-3Ar(x,£) dlGK{x,li) . 

Recall that the initial generating function was a third- 
order polynomial Gor(x) — giR x + 52R x'^ -\- (73R x^ with 
at least gsR < 0. Upon setting the right-hand side of 
Eq. p5|) to zero, we find the locally stable nontrivial 
fixed-point function 



Gr(x) = Ax 



(17) 



with arbitrary constant A, whereas the trivial solution 
G^ = is clearly unstable. This result is in fact valid 
to all orders in e [sj. Therefore, all running couplings 
5pr(^) ^ for p > 3, and we recover the results ((Tn|l . 
((TT|l . and ((T^ previously obtained from the reduced ac- 
tion @. As anticipated. A, the renormalized counterpart 
to (71, plays the role of a control parameter, albeit not a 
relevant one, since that would have to scale to infinity un- 
der renormalization rather than remain constant. We as- 
sume that A is a regular function of the initial rates; i.e., 
at fixed annihilation and coagulation rates we expand for 
(J ^ (Jc- A((t) ~ (fJc — cr) A'(ctc), since A(crc) = and 
A'(ctc) > 0. Indeed, the pure annihilation and coagula- 
tion model fixed points at cr = respectively correspond 
to A = e (since gi — 2^2) and A = e/2 {gi = (72) 0- 
Thus, in the PCPD inactive phase a < (Tc one should 
have A = 0(e) > as well. Notice that at cTc ^ 2fj, + ji' 
we have g2 = —3fj, — fj,' , hence this combination of an- 
nihilation and coagulation rates has apparently turned 
negative at the fixed point. One must therefore worry 
whether the physically accessible values of the reaction 
rates actually lie within the basin of attraction of the 
nontrivial fixed-point function ^^{x). 

Thus we now resume our discussion of the RG trajec- 
tories. On physical grounds one should expect that the 
flow would at least reach the line of flxed points encoded 
by Eq. (|17|l if giR is positive, since this corresponds to the 



inactive phase governed by the annihilation / coagulation 
fixed line. Yet for this to be true for any 51R > 0, the 
part of the scparatrix indicated by 11 in Fig. ^ must col- 
lapse onto the negative g2-axis, i.e., the separatrix should 
contain the invariant hypersurface gi = 0. This would in- 
deed lead to the standard RG flow picture: With initial 
values corresponding to the inactive phase the RG tra- 
jectories approach the annihilation / coagulation fixed 
line, critical initial conditions are defined by the separa- 
trix (which is unstable along one direction), and finally 
initial values corresponding to the active phase are to be 
found inside the unstable region with flow into the 'gully'. 

Yet it is easily seen that the zeros Xi with G-ii{xi,£) = 
are fixed by the partial differential equation (|15|) . and 
consequently also the intervals in which the function 
GR(a;, £) is respectively positive or negative. Therefore, if 
(i) giR > 0, there exists an open interval {xq — 0,xi > 0) 
wherein Gr{x,£) > 0, while Gb{xo_i,£) = 0. Since ini- 
tially (73R < 0, however, it follows that the fixed-point 
function (|17|) cannot be reached from the physically al- 
lowed region in parameter space. This is true for ar- 
bitrary values of A and for all —00 < a; < 00, even if 
5'2R > 0. This is rather astonishing because one would 
expect that at least deep in the inactive phase, where 
the production reactions can be neglected, Gn{x,£) 
Goo{x) — G^(x) as £ ^ 00 for all x. It is however ob- 
vious that we can assume this limiting relation to hold 
only in the first positive interval in an expansion of Gr (x) 
with respect to x. We conclude that the difference func- 
tion H{x) = G^(a;) — Goo (2;) must display an essential 
singularity at a; = 0, and its expansion in a series of x 
produces simply a zero. However, H{xi) — G^{xi). A 
qualitative discussion of the differential equation H15|l in- 
deed supports this assumption. Note here that GR(a;,.^) 
increases most significantly with £ at those values of x 
where its curvature is minimal, (ii) In the case giR < 
the initial function is negative for all x > 0. Hence the 
differential equation (|15|l does not provide any mecha- 
nism that could translate GR(a;, £) to positive values and 
finally to G^{x), at least in some x-intervals, without 
producing zeros of Gr(x) along the way. In summary, 
we find that our physical initial conditions at the critical 
point (with (73R < as well as g2R < for giR = A = 0) 
inevitably take the RG fiow from Eq. H15|) into the abso- 
lutely unstable region in Fig. Instead of reaching the 
stable fixed line (|17|l . we face runaway trajectories, to all 
orders in the perturbation expansion. 



IV. CRITICAL PROPERTIES AT THE 
UNPHYSICAL FIXED POINT 

Even though we have just seen that the critical fixed 
point gm, = A = and ^jr = e/2 is inaccessible to the 
RG flow trajectories starting at physical initial parame- 
ter values {i.e., positive reaction rates), let us neverthe- 
less explore the (hypothetical) ensuing critical behavior. 
Recall that in the PCPD, totally neglecting particle ex- 
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elusion, as encoded in the parameter Ai, suppresses the 
finite-density steady-state. Hence we retain this (appar- 
ently irrelevant) coupling, and moreover investigate how 
its RG flow towards zero becomes renormalized through 
fluctuations. From the explicit two-loop result for the 
associated renormalization constant Hll(l , we may imme- 
diately compute the anomalous dimension 



dlnZ 



7Ai 



Ai 



dn 



-3e + 3e2 In- + 0(e-'^) . (18) 



As is easily seen by investigating the RG equations for 
the particle density and its correlation function, al- 
ready completely determines the critical exponents here. 
This assertion is also confirmed by directly computing 
the renormalized equation of state (upon approaching the 
transition from the active side). For A ~ cr — ctc > 0, i.e. 
in the active phase, one finds that the steady-state den- 
sity vanishes as A ^ 0^ according to (0) ^ |A|'', with 



/3 = 



I + 7A1/2 



1 - 2e 



l + 31n-]e^ + Oie') , (19) 



while the two-point function correlation length (finite 
only in the active phase) diverges as ^ ~ where 



7Ai 

2 



1 + 



L^lnUoie^). (20) 



Precisely at the critical point (A 
decays asymptotically as {(j){t)) ^ 



0) the particle density 



t- 



with 



^=1 + 2^^1-3 3 4 

2 4 2 4 4 3 ^'^ ^ ' 

Since remarkably the anomalous fluctuation conrrec- 
tions to the critical exponents /3, v, and d are solely con- 
tained in the anomalous dimension (|18|l . we may elim- 
inate the latter to yield the following hyperscaling rela- 
tions, valid to all orders in e = 2 — d: 



25 + 1/v 



(3 = 25/{d-25) =dv-l 



(22) 



Here we have used the exact result for the dynamic ex- 
ponent z = 2 and the standard scaling relation f3 = zv 5 
(which also follows directly from the RG equations). At 
the critical dimension dc = 2, we infer the asymptotically 
exact scaling behavior 



A^O-: (^)^|A|(ln|A|)^ 

e~|Ari|ln|A|P 
A = 0: ((/)(t)) (1^^)3/2 _ 



(23) 
(24) 
(25) 



Aside from the fact that these exponent values are at 
odds with the presently available data from numerical 
simulations for the PGPD, they also lead to a serious 
contradiction, which confirms again that the fixed line 
(|17|l does not represent a physical system. First, we note 



that the positive value (72R = e/2 indicates the presence 
of particle aniicorrelations at this fixed point, precisely 
as in the pure binary annihilation or coagulation system. 
Next, recall that the equal-time density correlation func- 
tion of point-like particles consists of three contributions, 

(n(r,0n(r',i)) = (n(r,t))<5(r-r') 

+C{r,r';t) + {n{r,t)){n{r\t)) . (26) 

The first term here describes the particles' Poissonian 
self-correlations. C(r, r';t) represents the density cumu- 
lant, i.e., the connected correlation function of the den- 
sity fiuctuations, which is negative in the case of particle 
anticorrelations. Upon integrating Eq. (|26|) over the con- 
fining volume and dividing by the mean particle num- 
ber {N{t)), we obtain for a homogeneous state where 
C(r, r';i) = C(r — r';t) the following general expres- 
sion for the relative mean-square particle number fiuc- 
tuations: 



{N{t)) 



Cjt) 

m 



(27) 



with the mean density 4>{t) — {n{r,t)) and the spatial 
integral of the cumulant C{t). 

In the vicinity of a critical point these quantities scale 
as follows: 



0(t) - A t-f^'"" , C(t) - B t'^d^-m/-'^ 



(28) 



The amplitude A is of course positive, while _B < in 
the case of anticorrelations (such as in the PG univer- 
sality class), and i? > for positive particle correlations 
(as prevalent in the critical DP clusters). Gombining 
Eqs. (El and ^ yields 



(Nit)) ^'+a' 



(29) 



Thus, if di^ — /? > the second term dominates the right- 
hand side of Eq. H29|) asymptotically. For particle an- 
ticorrelations where B/A < 0, this would immediately 
contradict the positivity of the left-hand side. Conse- 
quently, the previously found scaling exponents which 
satisfy dv ~ (3 — 1 (exactly) are definitely inacceptable. 

This is in remarkable contrast to the results obtained 
for the pure pair annihilation or coagulation process 0, 
where the density and the integrated cumulant obey iden- 
tical asymptotic scaling behavior, 0(t) ^ C(t) ^ f'^^^. 
Hence no contradiction arises here provided \B\ < A. 
Neither are particle anticorrelations necessarily excluded 
at critical points: For even-offspring BAW that represent 
the PG universality class, a one-loop RG analysis at fixed 
dimension yields ly = 3/(10-3d) and /3 = 4/(10-3d) 0. 
Hence di^ — /3 < in dimensions d < d'^ = 4/3, which is 
precisely the borderline dimension (within the one-loop 
approximation) for the existence of the phase transition 
and the power-law inactive phase in this system. 
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V. CONCLUSIONS 

We have investigated the 'microscopic' field theory for 
the PCPD, as derived directly from the defining mas- 
ter equation, by means of the dynamic renormalization 
group. Although fluctuations generate an infinite chain 
of particle production processes A + A ^ {n + 2) A, the 
theory remains controlled and renormalizable in the in- 
active phase 0, where it is governed by the fixed point 
of the pure annihilation / coagulation model 0|- This 
is most elegantly seen by means of a functional RG ap- 
proach • In order to render the particle density finite 
in the active phase, we have incorporated site occupation 
restrictions following the methods developed in Ref. . 
On the mean-field level, this indeed leads to a continu- 
ous transition separating the active from the absorbing 
phase. However, a detailed analysis of the RG flow equa- 
tions shows that the action ((T)) does not adequately cap- 
ture the critical properties of the PCPD: (i) There is no 
stable RG fixed point in the physical region of parame- 
ter space , and instead one obtains runaway trajectories; 
(ii) the scaling exponents found at the RG fixed point in 
the unphysical regime violate the positivity of the mean- 
square particle number fluctuations. 

We emphasize that these statements in fact hold to all 
orders of the perturbation expansion and even apply to 
the nonperturbative 'exact' RG approach. This failure 
really resides in the starting field theory action itself, not 
its subsequent evaluation. Obviously, a crucial ingredi- 
ent was left out when the 'naive' continuum limit was 
performed. We can only speculate as to what the po- 
tential remedies might be, in part motivated by pictures 
from simulation studies, where positive particle correla- 
tions are observed both at the critical point and in the 
active phase: Perhaps one needs to explicitly introduce 
separate fields respectively for the positively correlated 
clustered particles and the solitary random walkers as 
coupled slow variables. The challenge then, however, is 
to write down a consistent coarse-grained field theory 
that correctly accounts for the internal stochastic noise 



generated by the reactions. To date, therefore, an apt 
field theory description of the PCPD remains an open 
and difficult problem. 

Finally, we remark that the same problems arise with 
the master equation field theory when the PCPD order 
parameter is coupled to a static, conserved background 
field js^. The above statements also apply to active- 
to-absorbing transitions in related higher-order reactions 
[TtI I2I 12^ For purely triplet reactions, 3^ mA 
(to = 0,1,2) and 3A {n + 3) A one readily derives 
a field theory representation corresponding to the action 
((Tj), essentially just raising the powers of (j) by one in 
the nonlinear vertices. In the inactive phase, the the- 
ory can again be analyzed by means of the functional 
RG analogous to Eq. (|15|l . leading to the upper critical 
dimension dc — 1, where the particle density decays ac- 
cording to {(j){t)) - {t-^lnt)^/"^ 0. At the hypothetical 
critical point one would obtain a slower critical decay 
{(/yit)) ~ t~^^^ (Inf)^/^. But once again, the critical fixed 
point cannot be reached by the RG flow starting at phys- 
ical parameter values. Moreover, the presence of a redun- 
dant parameter even questions the validity of identifying 
d = 1 as the critical dimension. Likewise, the expec- 
tation that in general fourth-order processes are merely 
described by mean-fleld scaling exponents may not be 
borne out by a correct treatment. 
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